Unconventional behavior of Dirac fermions in three-dimensional gauge theory 
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We study unconventional behavior of massless Dirac fermions due to interaction with a U(l) gauge 
field in three time-space dimensions. At zero chemical potential, the longitudinal and transverse 
components of gauge interaction are both long-ranged. There is no fermion velocity renormalization 
since the system respects Lorentz invariance. At finite chemical potential, the Lorentz invariance 
is explicitly broken by the finite Fermi surface. The longitudinal gauge interaction is statically 
screened and becomes unimportant, whereas the transverse gauge interaction remains long-ranged 
and leads to singular renormalization of fermion velocity. The anomalous dimension of fermion 
velocity is calculated by means of renormalization group method. We then examine the influence of 
singular velocity renormalization on several physical quantities, and show that they exhibit different 
behavior at zero and finite chemical potential. 

PACS numbers: ll.lO.Hi, ll.lO.Kk, 71.10.Hf 



I. INTRODUCTION 

Quantum electrodynamics of massless Dirac fermions 
defined in three time-space dimensions (QED3) has been 
extensively investigated for thirty years [1-23]. Different 
from its four-dimensional counterpart, QED3 is a super- 
renormalizable field theory, so its ultraviolet behavior can 
be well controlled. This gauge field theory is known to 
exhibit asymptotic freedom [10], which means the gauge 
interaction becomes stronger at lower energies. Due to 
this feature, intriguing non-perturbative phenomena are 
expected to occur in the low-energy regime. When Dirac 
fermions are strictly massless, the Lagrangian respects a 
continuous chiral symmetry. However, the strong gauge 
interaction may trigger fermion-antifermion vacuum con- 
densation, 7^ Oj which generates a finite fermion 
mass and induces dynamical chiral symmetry breaking 
(DCSB). Appelquist el al. first found that DCSB can 
occur when the fermion flavor N is smaller than some 
critical value Nc = 32 /tt^ Q. Motivated by this very 
interesting prediction, intense theoretical effort has been 
devoted to studying this problem [3-0, H-H^l- Despite 
some debate P,M 01; most of these analytical and nu- 
merical calculations agree that a critical value exists at 
roughly Nc w 3.5. A remarkable consequence of DCSB 
is that it leads to weak confinement Q . 

Apart from exhibiting many properties interesting in 
the context of elementary particle physics, QED3 also 
has extensive applications in condensed matter physics. 
Specifically, it serves as an effective low-energy theory 
of d-w&ve high-temperature superconductor [ij, [24| - [38j 
and some quantum spin liquid state [39|] . The happening 
of DCSB in QED3 corresponds to the formation of two- 
dimensional long-range antiferromagnetic order [l^ . [28l . 
nil US Recently, it is proposed that QED3 may 

be simulated on optical lattice [401, which provides an 
opportunity of measuring DCSB experimentally. 

When the fermion fiavor is beyond the threshold, N > 
Nc, no dynamical fermion mass can be generated. There 
is no fermion vacuum condensation, i.e., (iptp) = 0, and 



the continuous chiral symmetry is preserved. However, 
the absence of chiral condensation in the vacuum does not 
mean that the chiral symmetric phase is trivial. On the 
contrary, many highly nontrivial features can emerge in 
the symmetric phase of QED3. In particular, the gauge 
interaction is able to cause breakdown of Fermi liquid 
(FL). In 1973, Holstein et al. showed [HI that the un- 
screened transverse component of electromagnetic field 
in (3-f-l)-dimensional non-relativistic electron gas leads 
to unusual logarithmic specific heat, oc TlnT, which is 
apparently out of the scope of FL theory. This discov- 
ery have stimulated extensive investigations of non-FL 
behavior in various gauge theories, in the contexts of 
both condensed matter physics [13, l42| - [53j and particle 
physics [54,-58]. The non-FL behavior of massless Dirac 
fermions in QED3 has also been discussed [U HI, |3l| . In 
addition, QED3 can be used to describe some intriguing 
states, such as algebraic spin liquid state [29, - ,33., ,36] and 
algebraic charge liquid state 37]. 

In some many-particle systems described by QED3, 
there is a finite density of massless Dirac fermions. The 
finite fermion density is usually represented by chemical 
potential /i. An interesting question is how this chemical 
potential affects the physical properties of QED3. The 
impacts of finite chemical potential on DCSB have been 
addressed in [l3, [IE, HO] , where it is found that the crit- 
ical flavor Nc is lowered as chemical potential grows and 
DCSB is completely suppressed when chemical potential 
is sufficiently large. In the chiral symmetric phase with 
strictly massless Dirac fermions, the gauge interaction 
also leads to different properties at zero and finite chem- 
ical potential. For instance, the Dirac fermion dam ping 
rate behaves as oc ui^^^ at zero chemical potential [2l| 
and oc uj^^^ at finite chemical potential [22!]. 

In this paper, we consider chiral symmetric phase of 
QED3 and study unconventional behavior of massless 
Dirac fermions at finite chemical potential. We assume a 
relatively large fermion fiavor N so that the fermions are 
strictly massless. One important effect of a finite chem- 
ical potential is that it explicitly breaks the Lorentz in- 
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variance. As a consequence, the longitudinal component 
of gauge field develops an effective mass that is propor- 
tional to the chemical potential /z, which is analogous 
to the static screening of Coulomb interaction. However, 
the transverse component of gauge field remains massless 
due to the gauge invariance. Because of such difference 
in the longitudinal and transverse components of gauge 
field, the temporal and spatial components of fermion 
self-energy are no longer identical, which in turn gives 
rise to nontrivial renormalization of fermion velocity. 

We shall analyze fermion velocity renormalization by 
performing a renormalization group (RG) calculation 
|59l - l62l |. The fermion velocity remains a constant at 
zero chemical potential after including the gauge interac- 
tion corrections since the Lorentz invariance is preserved. 
However, it acquires strong momentum dependence after 
developing an anomalous dimension at finite chemical 
potential. The appearance of nonzero anomalous dimen- 
sion is a consequence of Lorentz symmetry breaking and 
gauge symmetry. Therefore, QED3 with a finite chemical 
potential is fundamentally different from that defined at 
zero chemical potential. We then evaluate specific heat, 
density of states (DOS), and compressibility of massless 
fermions using 7„ obtained in the RG analysis. After an- 
alytical and numerical calculations, we demonstrate that 
massless Dirac fermions exhibits unconventional, non-FL 
like, behavior at finite chemical potential. 

The rest of this paper is organized as follows. We define 
the Lagrangian and then perform RG analysis of fermion 
velocity renormalization in Sec. HIl A finite anomalous 
dimension of fermion velocity is obtained. We calculate 
specific heat in Sec. IHIl and DOS and compressibility in 
SecHVl In section |Vl we briefly summarize the results. 



the other hand, we consider the effective QED3 theory 
that is derived from a microscopic model of some con- 
densed matter system (for instance, Hubbard model of 
high-temperature superconductor) 2A^2^, the fermion 
velocity should be calculated from the band structure of 
the corresponding microscopic model. In the latter case, 
the fermion velocity vp is no longer equal to the velocity 
of light, c. Its magnitude is strongly material dependent, 
but is always much smaller than c |24l - l26| . In the present 
paper, we assume a constant velocity vp and study its 
renormalization due to gauge interaction at finite chem- 
ical potential. The main conclusion depends only on the 
momentum dependence of the renormalized velocity, but 
not on the concrete magnitude of the constant vp- 

In Euclidian space, the free propagator of massless 
Dirac fermions at zero ^ is 



°^ ' F • 



(2) 



At finite fx, it proves convenient to work in the Matsubara 
formalism and write the fermion propagator as 



Go(ia;„,k) 



{iujn + m)7o - vp^ ■ k ' 



(3) 



where the fermion frequency is iujn = i{2n + l)7r//3 — 
i{2n+ l)nT with n being integers. For notational conve- 
nience, we use fco to denote the imaginary frequency iujn, 
so that the fermion propagator can also be written as 

Go(fc) = r = 7o(fco+V;)+^i-7-k_ 



II. RENORMALIZATION GROUP ANALYSIS 
OF FERMI VELOCITY RENORMALIZATION 

The Lagrangian density for QED3 with iV-flavor Dirac 
fermions is given by 

^ - 1 

1=1 

where /i represents the chemical potential and vp is 
the constant fermion velocity. The Dirac fermions can 
be described by a four-component spinor field ip and 
■0 = T/jt-yO. The 4x47 matrices can be chosen as: 
7^ = ((T3,cri,(T2)'8'cr3, which satisfy the standard Clifford 
algebra {7^,71.} = 2g^,, with metric = diag(l, 1, 1). 
In this paper, we consider a large N and perform 
expansion. For convenience, we work in units with 
h — ks — I and restore them whenever necessary. 

It is now helpful to further remark on the physical 
meaning of fermion velocity vp. There are two ways to 
define the velocity vp. If we regard QED3 as a standard 
relativistic quantum field theory, the velocity of massless 
particles is simply the velocity of light, i.e., vp = c. If, on 



A. ^ = 

We first consider the case of zero fi. The key physical 
quantity is the fermion self-energy function, which can be 
used to calculate the RG flow of ifermion velocity vp [60I - 
l62j . Following the RG strategy presented in Ref . [g^ , we 
introduce two cutoffs Aq and Ai, with Ai being smaller 
than Aq. To the leading order of expansion, the 
one-loop fermion self-energy at zero fi is 

^^^^ ^ N A, (2^^-(fc3^^''^-('^) 

= loko^o + vp'jik,^!. (5) 

The full gauge field propagator D^u{q) is given by the 
Dyson equation 

D-^{q)^D^^y\q)+U,M), (6) 
with free gauge field propagator 

D^Kl) = ^,{9,.~'-f), (7) 
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in the Landau gauge. To the leading order oil/N expan- 
sion, the one-loop contribution to vacuum polarization 
tensor 11^^ is 



(27r)3 k^{q + kY 



(8) 



Here, it is convenient to define a = Ne'^, which is fixed as 
N taken to be large Q. Now the gauge field propagator 
becomes 



1 



g2 + n(g)V^"'' 



(9) 



where n(g) = ^ at zero temperature. 

According to the Dyson equation, G^^(fc) ~ Gq ^(fc) — 
the full fermion propagator is 



G(fc) 



1 



it - S(fc) 



1 



(1 - So) 70^0 + (1 - Si) VFjih 



(10) 



Apparently, (1 — Sq) correspond to the wave function 
renormalization, whereas (1 — Si) represents the prod- 
uct of the renormalization factors of wave function and 
fermion velocity. 

By inserting Eq. © into Eq. (O, it is easy to find that 
So = Si. Apparently, the temporal and spatial parts of 
fermion propagator are equally renormalized, which orig- 
inates from the Lorentz invariance of the system at zero 
/i. Therefore, the fermion velocity vp is a flow-invariant 
constant. 



B. fi^O 

At finite finite chemical potential, QED3 theory can 
exhibit qualitatively new features. First of all, there ap- 
pears a finite Fermi surface , which explicitly breaks the 
Lorentz invariance. In this case, the temporal and spatial 
components of fermion self-energy are no longer equiva- 
lent, so the fermion velocity may be singularly renormal- 
ized. 

To proceed, we first need to discuss the effects of a 
finite chemical potential fi on the effective gauge interac- 
tion function. These effects are reflected in the vacuum 
polarization function n(q, /x). It is technically hard to 
obtain an entirely analytical expression for n(g,/i), so 
we will derive an approximate Tl{q, fJ.) that captures the 
essential features of QED3 at finite /i. 

Under the Coulomb gauge condition, 9^11^ {q) — 0, the 
vacuum polarization tensor can be decomposed into two 
independent parts (63j : 



(11) 



where 



= S^oSa,^, (12) 
= 5p.(<5,, -^)(5,,, {*,j = l,2}. (13) 



They are orthogonal and related by 

A I R - A - Mil 



(14) 



Now, the full gauge field propagator Dfj^^{qo, q, /3) can be 
written as 



B„ 



q2 +n^(go,q,/3) g2 + nB(gro,q,/3)' 



(15) 

where the functions 11^ and 11 b are related to the tempo- 
ral and spatial components of vacuum polarization tensor 
n^i. by 



Ho 



He — n,« ^IIoo, 



(16) 
(17) 



with = 9o + q^- Based on these quantities, it is con- 
venient to write the fermion self-energy as 



S(fc) = Sa + Ss = 7o^oSo + VFlihT^i, (18) 



where 



5^0 — 5]^o + 5]so, 
Si = S^i + YiBi- 



(19) 
(20) 



In an earlier publication 22], the full analytical ex- 
pressions of polarization function n(g, ^) was derived. 
These analytical expressions are too complicated and can 
not be directly used. It is necessary to make proper ap- 
proximations. We first consider the temporal component 
noo(Q, A*). In order to simplify computations, we uti- 
lize the so-called instantaneous approximation [g^, i.e., 
Qo — 0, and write noo((?,/i) approximately as 



noo(q,At) = 



a|q| |q| 



(21) 



Apparently, the finite jj, serves as an energy scale: 
1100(9, m) behaves quite differently above and below 2/j,. 
In the low-energy (long wavelength) limit, q — 0, we 
have 



noo(q^O,//) 7^0. 



(22) 



As a consequence, the temporal part of gauge interac- 
tion becomes short-ranged due to static screening. In 
other words, the temporal component of gauge field ac- 
quires a finite effective mass that is proportional to fi. 
Such short-ranged interaction does not lead to any sin- 
gular contribution to fermion self-energy, and thus can 
be simply neglected. 
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We then consider the spatial component of polarization 
function 11^^(5, /i). Similar to its temporal counterpart, 
nii(9,/^) also exhibit different behavior above and below 
the energy scale 2/i. When |q| > 2/Lt, we still use the 
instantaneous approximation and have 



a\q\ 



(23) 



When |q| < 2^, Ilii{q,fi) vanishes at zero energy go = 0, 
so it is not appropriate to use the instantaneous approx- 
imation. The energy dependence of 11^^ (g, ^) should be 
explicitly maintained [22j. As pointed out in Ref. [23 |, 
the fermion self-energy is dominated by the low-energy 
region of go ^ |q| *C 2/j,. We notice this approximation 
is widely used in the calculations of fermion self-energy 
due to gauge interaction [41, 43, 44, 47, 48] and critical 
ordering fluctuation [g^]- In this region, Ilii{q, jj) can be 
approximated as 



n,;,(go,q,M) = i^n- 

Zn |q| 

At the lowest energy qq — 0, we have 
^ii{qo = 0,q,M) = 0. 



(24) 



(25) 



This fact implies that, the transverse component of gauge 
interaction remains long-ranged even when the dynam- 
ical screening effect is taken into account. It also ex- 
plains why the instantaneous approximation can not be 
used in this region. Physically, the long-range property 
of the transverse gauge interaction is protected by the 
local gauge invariance. 

The fermion self-energy is given by 



^« d^q (^ 



TV A, (2^)3 {k - q)^'^-q^ 



d^q -7o(fco-go) + Wi=^7-(k-q) 



NJk, (2^)3 



(fc - qf 



n,; 



^° d^q -2vFlikiql - 2vFl2k2ql 1 



NJa, i2nf 



q^(fc — g)^ 



(26) 
(27) 



= (7oA;oSbo + W_F7ifciS-Bl)fir 

-|-(7ofcoEBo + VF^tkiY,Bi)s 

According to the general RG scheme [s^, one needs to 
integrate out the high-energy degrees of freedom (fast 
modes) step by step, until eventually reaching the lowest 
energy. Since 11^^ behaves differently at high and low 
energies, the self-energy should be calculated separately. 

For small external momentum fc <C Ai, we are allowed 
to make the approximation [60J, (fc — g)^ g^, and write 
the first term of Eq. (|27)) as 

r-Ao 



(S_B1 — EBo)fir 



d^q 1 



1 



2a 

^A, WrqHq^+I^,^) 
2a d^q 1 1 



(28) 



From RG theory, we know that the fermion velocity 
can receive singular renormalization only when the self- 
energy contains a logarithmic term. One can check that 
there is no such term in the regime where <C g^. 
Nevertheless, a logarithmic term emerges as one goes to 
the low-energy regime where ^ g^ . After neglecting 
g^, it is straightforward to obtain 



Bl 



S_Bo)fir 



2a d^q 1 1 

A^Ja^ (2^7 2M 

ttN Ai 



(29) 



Paralleling the above analysis, the second term of Eq. 
can be computed in an analogous manner: 



-2^ 

N 


f cPq 


qh 1 


J (2^)3 


q2(fc-g)2g2 + n,. 


-2^ 


f d^q 


qh 1 


N 


J (2^)3 


q2g2 n,, 


-2^ 
N 


f d^q 


qh 8 


J (2^)3 


q^g-^ a q 


2 


Ai 





(30) 



Now the total self-energy is 



(Ssi — T,Bo)fir + C^Bl — ^Bo)s 

4 , Ao 

^^"aT 



— In^ 



2 . Ao 
iV^^^AT- 



(31) 



This result is valid only when Aq > Ai > 2//. As the 
energy scale decreases below 2fi, Uu should be replaced 
by its low-energy expression. In this case, the first term 
of fermion self-energy becomes 



Bl 



SBo)fir 



rf3g 1 



1 



Nj^^ (27r)3g2g2 + n 
2a f^" d^q 1 1 



TV A, i^^rq'q' + ^^ 



(32) 



In order to simplify calculations, we can divide momenta, 
frequency, and chemical potential by a to make all these 
variables dimensionless. However, for notational conve- 
nience, we still denote Ao,i/q; as Ao,i. Now, the above 
equation is rewritten as 



Bl 



SBo)fir — -77 



^° d^q 1 



1 



TV A, (27r)3g2,2 + Ji.^ 



(33) 



In the low-energy region go <C |q| ^ 2/i, this integral 
becomes 



(S 



Bl 



SBo)fir 



2 

TV 



dqo 



d\q\ 



(277)2 |q|3 + 



5 



2 f f d\q\ 1 



27r)2 |q|3 + ^q. 



X0(A2-q2)0(q2_A?) 

2_ f^' /-^" rf|q| 

ivio '^°A, (2^)2|qP + ^9o■ 



l 



(34) 

It seems easier to first integrate over qo and then integrate 
over q, 



(S_Bi — E_Bo)fir 



4n 



In 1 



„ ^ In ( ^ ) + others terms. 
Stt^A' V Ai 



(35) 



Analogously, we have 

(S_B1 — SBo)sec = - 



—4— 7 In 1^1+ others terms. 
Stt'^N V Ai 



(36) 



We finally obtain the total contribution 



Here we keep only the logarithmic term, which survives 
at the lowest energy and corresponds to the stable fixed 
point produced by the long-ranged transverse gauge in- 
teraction. 

Summarizing the above results, we obtain the following 
RG equation 



where 



1 



71- 



(38) 



(39) 



87r2iV' 

This equation has the following solution, 

^)J.cxF^ (40) 

where the constant 7„ defines an anomalous dimension 
of velocity vp and the dimensionless momenta k actually 
corresponds to fc/A with A being the UV cutoff. Ob- 
viously, the constant fermion velocity becomes strongly 
momentum dependent due to the long-range transverse 
gauge interaction. The fc-dependence of fermion velocity 
is shown in Fig. ([T]) at both zero and finite fj,. If we define 
vf{Aq) as the bare velocity, which is taken to be unity 
in the figure, then the effective velocity vp{k) decreases 
monotonically when k is lowering. It eventually vanishes 
as fc — J> 0. 



FIG. 1: Momenta k is rescaledr by ultraviolet cutoff, k — >■ 
fc/A. The upper curve represents the constant velocity at jj. — 
0. The lower curve shows the fc-dependence of renormalized 
velocity at /J. 0. 



After this renormalization, the energy spectrum of 
Dirac fermions becomes 



e{k) oc fci+''" 



(41) 



This type of unusual velocity renormalization is a char- 
acteristic feature of Dirac fermion systems, and has 
been studied extensively in the contexts of rf-wave high- 
temperature superconductor [l^, HH, [G^I and graphene 
(goI l65l . |66| . The finite anomalous dimension of fermion 
velocity distinguishes the system at finite /i from that at 
zero /i. An interesting issue is how this anomalous dimen- 
sion affects the observable quantities of Dirac fermions, 
which will be addressed in the following two sections. 

It is interesting to note that the anomalous dimension 
7i, is a universal constant for any given flavor N. Al- 
though jy is finite only at finite /z, it does not explicitly 
depend on /i. It is the finiteness, rather than the precise 
value, of /i that is really important. A finite 7„ is gen- 
erated as long as /i becomes finite, no matter how small 
it is. Therefore, the ground state changes fundamentally 
once one moves away from the neutral Dirac point. 



III. SPECIFIC HEAT OF DIRAC FERMIONS 

According to Landau FL theory, a weakly interacting 
fermion system is in one-to-one correspondence with a 
non-interacting fermion gas. There is a sharp Fermi sur- 
face in a FL with well-defined quasiparticles existing in 
the low energy regime. The properties of a FL can be 
manifested in a variety of physical quantities, such as 
spectral function, specific heat, DOS, and susceptibility. 
The Coulomb interaction does not destroy the stability of 
FL state in normal metals, because it is always statically 
screened by the collective particle-hole excitations. How- 
ever, non-FL behavior would emerge when there is some 
kind of long-range gauge interaction. We now would like 
to examine the corrections of gauge interaction to several 
physical quantities of massless Dirac fermions. We will 
consider specific heat in this section, and consider DOS 
and compressibility in the next section. 
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In order to calculate specific heat, we first need to 
calculate the free energy. In the following, to facili- 
tate calculations we make the rescaling transformations: 
T, — ;> T/A, fc/A, /i/A, where A is ultraviolet cutoff. 
When n — 0, the specific heat due to gauge interaction 
has already been studied in previous works [s^jUHl, which 
show that Cv oc T^. 

The free energy density will be computed using the 
methods given in Ref. i6T]. The partition function is 



^ - n y'2?H<n(k)]2?[V'o,n(k)]e^, 

n,k,a 

where the action is expressed as 

S = ^[-*Vln(k)]i?a,p[V'p,n{k)], 



(42) 



Ti.k 



D 



P[-{^n - - ivpl^l ■ k] 



(43) 



Employing the functional integral formula for Grassmann 
variables, 



(44) 



P[?7^]2?Me'' = detZ? 
it is easy to obtain 

InZ = ^lndet£)(n,k) 

Ti,k 
n,k 



n.k 



vl{kl + 



n.k 



+ \u[p\u^l + {e{k)+^lf)\} 
Applying the identities. 



(45) 



1 e^ + {2n+l)^n 

= ln[(2ji + l)27r2+/32(aj±^)2]^ 

1 7r[cot(7ra::) — cot(7ri/)] 



E 



n— — oo 



{n - x){n - y) 
we are left with 



+ ln[l + (2n + ifiT^] 
(46) 
(47) 



InZ = 2V 



+ lnfl + e-'3(^W+'') 



/3e(fc)+ln(l + e-'3(^W-^)) 



(48) 



after some straightforward algebra. Since the free energy 
density f — ^ — —-^y In Z, we have 



/(r,7.,M) = -2ry ^In 



1 + exp 



j(fc) ±^ 
T 



(49) 



where the zero-point energy (the first term in InZ) has 
been discarded. The fermion specific heat can be ob- 
tained from free energy, namely 



-T 



QT2- 



(50) 



In the absence of gauge interaction, the fermion ve- 
locity takes its bare value and the anomalous dimension 
7i, = 0. Therefore, the free energy density is simply 



-2T 



47r2 



In 1 



^)+ln(: 



It is easy to integrate over k, and get 



3 r 



Li.-? 



Lis -e T 



(51) 



(52) 



where ^[^{z) is polylogarithmic function. The specific 
heat has the following form. 



-4/iT 



1 + e^T 

L 

Li2 ( -e^T 



1 + e^-' 



13 (-e^) +Li3 (-e 



Li2 



}.(53) 



It is obvious that the specific heat = 2CMt'^ qc 
at /i = 0. At finite /z, the transverse gauge interac- 
tion induces an anomalous dimension for fermion veloc- 
ity, which modifies the Dirac fermion energy spectrum. 
Now the free energy density becomes 



/(Ai,7i.) 



-2T /^In 
27r 



(l-l-e T ja + e T \ 



(54) 



where ?7 = 1 + 7t,. It is convenient to define x — , and 
write the free energy as 



f{f^,v) = / dxx 'J In 

TTI] 



/../-ln[(l + e-J^)(l + e- 



(55) 



Making derivatives of /(^, rj) with respect to T, we obtain 



Cv = 



dxx' 



{x + fj,)'^e~ 



1 I ^ 

i + e T 



(a; — /i)^e 



(56) 

which is complicated and will be evaluated numerically. 

Note the UV cutoff A does not qualitatively affect our 
basic conclusion, which allows us to set A = 1. The 
specific heat Cy explicitly depends on both chemical po- 
tential fi and temperature T. Its T-dependence shown 
in Fig. ([2]) for iV = 4. When chemical potential fi = 0, 
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FIG. 2: 7u = g^ijv • The curves from top to bottom are 
Cv(/i = 10"^), Cvifi = 10"'^), Cv(m = 10"*), and = 
2^T^, respectively. 



The curves from top to bottom are 



FIG. 4: 7„ 

Cv(m = 10-''), Cvin = lO"*^), Cv(/i = 10"*), and 
^^^T^ respectively. 




FIG. 3: 



7« 



57„. The curves from top to bottom are Cv{l^ 



10--*), Cvifi 
respectively. 



10-"), Cv(p = 10-**), and C'i 



9C(3) y2 



the specific heat behaves as Cy oc T^. At finite /i, the 
specific heat deviates from the curve. The deviation 
becomes more significant for larger /i. At lovif tempera- 
ture, the specific heat can be approximately written as 
power- law, Cy oc , where the exponent 5 is a function 
of fi. Such unconventional non-FL behavior arises from 
the anomalous dimension of fermion velocity, 7^,, which 
is generated by the long-ranged transverse gauge interac- 
tion at finite fi. In particular, the deviation occurs once 
fi becomes finite. This implies that the ground states of 
QED3 are very different at zero and finite fi. 

In many works on non-FL behavior caused by singular 
interactions, the specific heat is expressed in a logarith- 
mic function [Hill E H HI HI , ie., Cy oc TlnT. 
Note that this logarithmic expression does not contra- 
dict our results. Actually, the power-law specific heat 
presented here amounts to a summation of all powers of 
InT [5l,[53,[63,[6i. 



In the above calculations, the anomalous dimension 7^ 
plays an essential role. However, it is obtained by adopt- 
ing certain approximations. To examine the reliability of 
our results, we suppose 7„ = g^Jipf becomes 57i, or 7i,/5 
after including higher order corrections, and show the 
corresponding results in Fig. ^ and Fig. ([4]). Appar- 
ently, the basic conclusions are independent of the precise 
value of anomalous dimension. 



IV. DENSITY OF STATES AND 
COMPRESSIBILITY 

We now turn to the interaction corrections to DOS and 
compressibility. At /i = 0, the propagator of massless 
Dirac fermion is 



1 



(57) 



After analytical continuation, iuj u + iS, we have the 
following retarded propagator 



1^70 + VFJ ■ k 

j'^ — Vpk'^ + isgn{uj)6' 



(58) 



The corresponding spectral function is given by 
A(w,k) = -ilmG"'(w,k) 

TT 

= sgn(w)(a;7o - vp-f ■ k)S{uj^ ~ wf^k^), (59) 
which then gives rise to the DOS, 

= nJ -^Tl'{7„ImG'"(u.,k)} 

^ ' (60) 

Apparently, the DOS vanishes at the Fermi level, a; = 0. 
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At finite chemical potential /z, the fermion propagator 
becomes 



which leads to 



G(ia;,k) 



{iuj + fi)jQ + vpj ■ k 
{iuj + — Vpli^ 



(61) 



The DOS can be calculated similarly, with the expression 



p{uj) = N 



(2^)2 
N{uj + ii) 



Tr{7oImG"'(tJ,k)} 



(62) 



which approaches a constant proportional to /i as a; — > 0. 
We are interested in the gauge interaction corrections to 
the above expressions. At finite /i, the fermion velocity 
vp becomes A:-dependent, v'p oc fc'''". Using this expres- 
sion, we find that, 



AN'k{uj + fi) 
4NTT{bJ + ii) 

N{uj + ii)^ 



dk' 

(2^ 



(63) 



The effect of gauge interaction is reflected in the nontriv- 
ial exponent of DOS. 

Since the difference between zero and finite chemical 
potential is of primary interest in our work, we only 
consider zero-temperature compressibility, which is an 
important quantity describing the electronic properties 
of any interacting system. In its original meaning, the 
compressibility is defined as k = dV/dF with V 
and F being the volume and compressing force, respec- 
tively. However, in practical many-particle calculations, 
it is more convenient to define the compressibility as 
K = dn/dn, where n is the number of particles per area 
[ggI [70| . The compressibility of massless Dirac fermions 
vanishes at zero chemical potential, fi — 0. This behav- 
ior will be changed at finite density, i.e., /x ^ [66l[70j. 
In order to obtain n, we only need to calculate the fi- 
dependence of particle number n. 

Using the DOS expressed in Eq. (|62|) . the particle 
number in the absence of gauge interaction is 



dujp{uj) = 



" A.(. + ,)^_A^^ (64) 



2w|,7r' 



dn N p. 

dpL VpTT 



(65) 



After including the gauge interaction, the DOS is given 
by Eq. (|63l) . Now the particle number becomes 



n — diopiio) — I y 

which then yields 



duj — 



2vj,Tr 



K = 



(66) 



(67) 



Once again, the effect of gauge interaction is reflected in 
the exponent. 



V. SUMMARY 

In summary, we have studied the effects of a finite 
chemical potential on the behavior of Dirac fermions in 
QED3. At zero chemical potential, there is no fermion 
velocity renormalization. At finite chemical potential, 
the longitudinal gauge interaction becomes short-ranged, 
but the transverse gauge interaction remains long-ranged 
and leads to singular velocity renormalization. An ex- 
plicit calculation shows that a finite anomalous dimen- 
sion of velocity is generated and gives rise to unconven- 
tional properties in some physical quantities, including 
specific heat, DOS, and compressibility. Therefore, the 
massless Dirac fermions behave quite differently at finite 
and zero chemical potential. This difference, together 
with the difference in fermion damping rate |2ll |22| , in- 
dicates that the ground state of QED3 is fundamentally 
changed once the chemical potential becomes finite. 
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